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NUCLEAR IMAGING SYSTEM USING SCINTILLATION BAR DETECTORS AND 
METHOD FOR EVENT POSITION CALCULATION USING THE SAME 

5 BACKGROUND OF THE INVENTION 

1. Field Of The Invention 

The present invention generally relates to nuclear 
medicine, and systems for obtaining nuclear medicine images 
of a patient's body organs of interest. In particular, the 
10 present invention relates to a novel detector configuration 
for nuclear medical imaging systems that are capable of 
performing either positron emission tomography (PET) or 
planar and single photon emission computed tomography 
(SPECT) . 

15 

2. Description Of The Background Art 

Nuclear medicine is a unique medical specialty wherein 
radiation is used to acquire images that show the function 
and anatomy of organs, bones or tissues of the body. 

20 Radiopharmaceuticals are introduced into the body, either by 
injection or ingestion, and are attracted to specific 
organs, bones or tissues of interest. Such 
radiopharmaceuticals produce gamma photon emissions that 
emanate from the body. One or more detectors are used to 

25 detect the emitted gamma photons, and the information 

collected from the detector (s) is processed to calculate the 
position of origin of the emitted photon from the source 
( i.e. , the body organ or tissue under study) . The 
accumulation of a large number of emitted gamma positions 

30 allows an image of the organ or tissue under study to be 
displayed. 

Two basic types of imaging techniques are PET or 
"coincidence" imaging, and single photon imaging, also known 
as planar or SPECT imaging. PET imaging is fundamentally 
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different from single photon imaging. In PET, events are 
detected from the decay or annihilation of a positron. When 
a positron is annihilated within a subject, two SllKeV gamma 
rays are simultaneously produced which travel in 
5 approximately opposite ( i.e. , 180°) directions. Two 

scintillation detectors are positioned on opposite sides of 
the patient such that each detector will produce an 
electrical pulse in response to the interaction of the 
respective gamma rays with a scintillation crystal. In order 

10 to distinguish the detected positron annihilation events 

from background radiation or random events, the events must 
be coincident ( i.e. , both occur within a narrow time window) 
in each detector in order to be counted as "true" events. 
When a true event is detected, the line connecting the 

15 positions of the two points of detection is assumed to pass 
through the point of annihilation of the positron within the 
subject being imaged. 

By contrast, single photon imaging, either planar or 
SPECT, relies on the use of a collimator placed in front of 

20 a scintillation crystal or solid state detector, to allow 

only gamma rays aligned with the holes of the collimator to 
pass through to the detector, thus inferring the line on 
which the gamma emission is assumed to have occurred. Both 
PET and single photon imaging techniques require gamma ray 

25 detectors that calculate and store both the position of the 
detected gamma ray and its energy. 

Present day single photon imaging systems all use large 
area scintillation detectors (on the order of 2000 cm 2 ) . 
Such detectors are made either of sodium iodide crystals 

30 doped with thallium (Nal(Tl)), or cesium iodide (Csl). 

Scintillations within the Nal crystal caused by absorption 
of a gamma photon within the crystal, result in the emission 
of a number of light photons from the crystal. The 
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scintillations are detected by an array of photomultiplier 
tubes (PMTs) in close optical coupling to the crystal 
surface. Energy information is obtained by summing the 
signals from the PMTs that detected scintillation photons, 
5 and position information is obtained by applying a 

positioning algorithm to the quantitative signals produced 
by the PMT array. The original gamma- ray camera is 
described in U.S. Patent No. 3,011,057 issued to Hal Anger 
in 1961. 

10 The Csl camera is typically used with either a single 

silicon-based photodiode detector or an array of silicon- 
based photodiode detectors, which detect scintillation 
events emitted from the Csl crystal. Csl crystals are used 
where the relatively low cost, ruggedness and spectral 

is response of the Csl crystal are desired in favor of 
alternative crystal materials such as Nal . 

In PET imaging, scintillation crystals with short 
response times are required in order to properly detect the 
coincidence events with high timing resolution. Typically 

20 such crystals are chosen from among materials such as Nal, 
BGO, LSO and BaF 2 . Detectors coupled to such crystals can 
be an array of PMTs, a single "position-sensitive" PMT ("PS- 
PMT" ) , or fast-response silicon-based photodiodes such as 
avalanche photodiodes . 

25 Because the conventional Anger camera uses a thin 

planar sheet or disk of scintillation crystal material, it 
is necessary to cover the entire field of view of the 
crystal with light detectors such as PMTs or photodiodes. 
Additionally, the sampling capability of such scintillation 

3 0 crystals could be improved by increasing the number of gamma 
photons emanating from an imaging subject that are absorbed 
by the crystal, and consequently increasing the number of 
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scintillation events that can be detected for use in 
constructing an image. 

The bar detector is a specific configuration of 
scintillation detector that has been used in astronomical 
5 and high energy physics applications. The bar detector 

consists of an elongated scintillation crystal bar having a 
relatively small cross section. A photosensor such as a PMT 
is optically coupled to each end of the bar. The light from 
a gamma photon event within the scintillation crystal volume 

10 is detected by the two PMTs . The timing or signal 

information can be used to determine the location of the 
event in the bar. Additional bars can be placed next to 
each other for two-dimensional detection. 

An example of a proposed design for a PET detector 

15 module using a bar detector is given in Moses et al . , 
"Design Studies for a PET Detector Module Using a PIN 
Photodiode to Measure Depth of Interaction," IEEE 
Transactions on Nuclear Science NS-41, pp. 1441-1445 (1994) , 
incorporated herein by reference in its entirety. According 

20 to this design, a scintillation bar is coupled at one end to 
a PMT, and at the other end to a photodiode, in order to 
measure the depth of interaction (DOI) of the scintillation 
event in the bar . 

In past bar detector experiments for physics and 

25 astronomy, Nal (and sometimes Csl) bars of up to 100 cm were 
used to detect gamma photons of up to 10 MeV. Positional 
resolution within the bar ranged from 1.5 cm at 200 keV to 2 
cm at 10 MeV, although worse resolutions were reported. An 
energy resolution of 9.4% and a timing resolution of 10 ns 

30 at 662 keV and a 100 cm Nal bar were reported by a physics 
group for a balloon borne gamma telescope. Energy 
resolutions from other experiments were higher for the same 
energy gamma photon. These studies have cited geometry, bar 
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size, light attenuation coefficient and electronic noise as 
the major factors in determining the spatial and energy- 
resolution of bar detectors. However, the performance of 
bar detectors as designed in the prior art is insufficient 
5 for use in medical imaging applications. 

It would be desirable to use a bar detector in medical 
imaging applications. First, the bar detector can achieve 
the same field of view as a conventional Anger camera with a 
significantly smaller number of PMTs . This is because a 

10 scintillation event occurring anywhere within the length of 
the bar can be detected by the light sensors optically- 
coupled to the ends of the bars using a DOI calculation. 
Thus, photodetectors do not need to be placed over the 
entire field of view as in the conventional Anger camera. 

15 Additionally, multiple bars can be placed next to each other 
in a two dimensional array, creating a detector module than 
can localize the gamma photon interaction in a three 
dimensional volume, allowing for depth of interaction (DOI) 
decoding capability. This DOI information can enhance the 

20 spatial resolution of a positron emission tomograph system 
by improving sampling characteristics of the detector 
system, as the total thickness of the module can be larger 
than the thickness of a single crystal and thus can interact 
with a higher percentage of gamma photons emanating from a 

25 subject. Additionally, for SPECT systems, the bar detector 
module presents a very economical modality for constructing 
a camera with substantial advantages over a conventional 
planar camera (see Figs. 3a-3b) , including high throughput 
and elimination of the necessity for gantry motion. 

30 

SUMMARY OF THE INVENTION 

The present invention solves the existing need by 
providing a gamma camera having a scintillation detector 
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formed of multiple bar detector modules. The bar detector 
modules in turn are formed of multiple scintillation crystal 
bars, each being designed to have characteristics to achieve 
a sufficient spatial resolution for nuclear medical imaging 
5 applications. 

According to another aspect of the invention, the bar 
detector modules are arranged in a three-dimensional array, 
where each module is made up of a two-dimensional array of 
bar detectors with at least one photosensor optically 
10 coupled to each end of the module. Such a camera can be 
used for both PET (coincidence) and single photon imaging 
applications . 

According to another aspect of the invention, a bar 
detector gamma camera is provided, which utilizes an 

15 improved positioning algorithm that greatly enhances spatial 
resolution in the z-axis direction ( i.e. , the direction 
along the length of the scintillation crystal bar. 

In particular, the present invention provides a gamma 
camera, including a number of bar detector modules arranged 

20 in a three-dimensional array, each of the modules including 
a number of elongated scintillation crystal bars, each bar 
having two end surfaces, the number of bars being arranged 
in a two-dimensional array with respect to the end surfaces, 
and at least two photosensors, each optically coupled to a 

25 respective end of the module. 

According to another aspect of the invention, a gamma 
camera is provided, including a number of bar detector 
modules, each including a number of elongated scintillation 
crystal bars, each bar having two end surfaces, the 

30 plurality of bars being arranged in a two-dimensional array 
with respect to the end surfaces, and at least two 
photosensors, each optically coupled to a respective end of 
the module, for detecting a scintillation interaction of a 
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gamma photon with one of the bars; and a position calculator 
for determining the spatial location of a detected 
scintillation interaction in the elongated dimension of a 
scintillation crystal bar, according to the formula: 



z = arg min 



where z is the elongated dimension of said bar, 



10 R = 



E l+ E 2 



Ei = the total energy detected at a first end of said bar, 
E 2 = the total energy detected at a second end of said 

bar, 

li R = the mean of ratio R at a given location z, and 

q r = the variance of the ratio R at a given location z. 

20 BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will become more clearly understood from 
the following detailed description in connection with the 
accompanying drawings, in which: 

FIGs . 1(a) -1(c) are isometric views of various 
25 configurations of scintillation bar detector modules for use 
as building blocks in a gamma camera, according to the 
present invention; 

FIGs. 2 (a) - 2(d) are graphs illustrating point spread 
function profiles for test points in scintillation crystal 
30 bars constructed according to the concepts of the present 
invention; 

FIGs. 3(a) - 3(c) are perspective views of various 
three-dimensional configurations of bar detector modules for 
use in gamma cameras according to the present invention; 
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FIG. 4 is a graph showing spatial resolution as a 
function of scintillation light yield for a bar detector 
constructed according to the present invention; 

FIG. 5 is a circuit block diagram of a position 
5 calculation circuit for calculating spatial location of 
scintillation events in a bar detector gamma camera 
according to the present invention, using position- sensitive 
photomultiplier tubes (PS-PMTs) as photosensor elements; and 

FIGs. 6(a) - 6(b) are graphs illustrating the effects 
10 of various surface treatments of scintillation crystal bars 
according to the present invention, with respect to spatial 
and energy resolution. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

15 Referring to Figs. 1(a) -1(c), according to one 

preferred embodiment of the invention, a bar detector module 
is constructed by bundling together a number of 
scintillation crystal bars in a two-dimensional array in the 
X-Y direction as shown. A nuclear medical imaging gamma 

20 camera then can be constructed by aggregating a number of 
modules together in a preferred geometry to establish a 
desired field of view for medical imaging applications. In 
the examples illustrated, the bars are arranged in a 7x7 
square; however, many other configurations are possible and 

25 contemplated according to the invention. The bar detector 
module is optically coupled at each end to at least one 
photosensor. As shown in Fig. 1(a), each module may be 
optically coupled to a position-sensitive photomultiplier 
tube (PS-PMT) . The PS-PMT is a known photosensor that has 

30 the ability to identify the X-Y location of light photons 
incident on its detection surface and produce an electric 
signal indicative of such X-Y location, in addition to 
indicating the intensity of the incident light photons. 
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Similarly, as shown in Fig. 1(c), the bar detector module 
may be optically coupled to a pair of photodiode arrays, 
where each photodiode in the array corresponds to an 
individual scintillation crystal bar in the module. In each 
5 of these arrangements, the signal intensity information from 
the photosensors at each end of the module is used to 
calculate the DOI information, i.e., location or position 
information along the Z-axis. 

Alternatively, as shown in Fig. 1(b), an array of PMTs 

10 may be optically coupled to the ends of a number of bar 
detector modules aggregated together in a particular 
geometry. The PMT array is similar in function to the 
conventional PMT array provided with a single scintillation 
crystal sheet. According to the invention, however, the 

15 position information provided by both PMT arrays is used to 
calculate DOI position information in the Z direction, as 
well as calculating X-Y position information. 

The X and Y directional spatial resolution is mainly 
determined by the cross -sectional area of the bar. However, 

20 there are many factors affecting the performance of the Z 

directional resolution such as cross-sectional area, surface 
treatment of the bar surfaces, and light output and length 
of the bar. A key factor among these is the length of the 
bar. The effect of bar length on Z-axis spatial resolution 

25 is shown in Figs. 2 (a) -2 (d) . In these graphs, the cross- 
section of the bar was simulated to be 4x4 mm 2 , with 
respective lengths of 5 cm, 10cm and 20cm for Figs. 2(a) , (c) 
and (d) , with the long surfaces being ground. A statistics- 
based positioning algorithm, described in detail below, was 

30 employed except in case (b) , where a conventional centroid 
positioning algorithm was used. In all cases, the 
photosensor quantum efficiency and scintillation photons per 
gamma event were set to 0.25 and 2 0,000, respectively. 
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The mean full width half maximum (FWHM) of the point 
spread function of test points along the z-direction was 0.9 
mm, 1.7 mm and 3.8 mm for the 5cm, 10cm and 20cm length 
bars, respectively, as shown in Figs. 2(a), 2(c) and 2(d), 
5 using the statistics-based positioning algorithm described 
below. In Fig. 2(b), a 5 cm long bar was used with a 
conventional centroid positioning scheme. FWHM is a 
technical term used to describe a measurement of the width 
of an object in a picture, where that object does not have 

10 sharp edges, as is the case with light photons spreading out 
from a scintillation event in all directions. The light 
distribution generally has a Gaussian profile. Because the 
"width" of the profile extends to infinity, it is not 
possible to describe the extent of such profile. An 

15 alternative way to describe the approximate size of the 

object is to specify the width across the profile when it 
drops to half of its peak, or maximum, value. This 
corresponds to the full width of the object image at half 
maximum value, or FWHM. Thus, smaller values of FWHM 

20 correspond to higher spatial resolutions. Fig. 2(a) shows 
that for a 5 cm long bar using the novel statistics-based 
positioning algorithm described below, a very good FWHM 
value of 0.9 mm was achieved. 

By selecting the orientation, length and the size of 

25 the cross-sectional area of the scintillation bar, diverse 
imaging systems can be implemented with desired system 
configuration, performance and cost. Figs 3(a) -3(d) show 
examples of such applications. For instance, the 
configuration in figure 3 (a) provides a low-end but quite 

30 economical system by using less photosensors per field of 
view (FOV) area, while the configurations shown in Figs. 
3(b) and 3(c) can be used for high resolution systems. 
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The light yield property of a scintillator detector 
also is important in determining the spatial and energy 
resolution of any proposed system. Prior art bar detectors 
exhibited insufficient spatial resolution (about 1cm FWHM) 
5 for nuclear medicine imaging due to the poor light yield 
properties of the scintillation crystal. 

Fig. 4 shows Monte-Carlo simulation results of spatial 
resolution as a function of the number of scintillation 
photons per gamma event . The mean FWHM value was calculated 

10 from a number of test points along a 20 cm long by 4x4 mm 2 
cross-section bar detector as a function of scintillation 
photons. The quantum efficiency (QE) of the photosensor was 
assumed to be 0.25. The results conform to the theoretic 
fact that spatial resolution is inversely proportional to 

15 the square root of the total light output. 

In order to achieve sufficient spatial resolution for 
nuclear medicine imaging applications, i.e. , at least -4-5 
mm FWHM, the light yield of scintillation crystal must be 
higher than about 12,000 photons per event, as shown by the 

20 results depicted in Fig. 4. This is equivalent to 

scintillating an Nal crystal with -3 00 keV gamma photons and 
collecting the scintillation light photons with a typical 
PMT array having 0.25 QE . Therefore, high light yield 
scintillation materials such as Nal, LaCl3, LaBr3 or LSO are 

25 preferred for high resolution medical imaging applications. 

LaCl3 and LaBr3 are newly developed crystals with very high 
light output and extremely fast decay times. By utilizing 
these new scintillator materials, the performance of the bar 
detector can be greatly improved over prior art systems, 

30 allowing for its use in nuclear medicine imaging. Physical 
characteristics of the scintillator materials are listed in 
table 1. 
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Table 1. Physical properties of inorganic scintillators with high light 
output . 





Nal 


LaC13 


LaBr3 


LSO 


Density (g/cm J ) 


3 .67 


3 . 86 


5.29 


7.4 


Z e ff 


51 


50 


47 


66 


Light output 
(ph/MeV) 


39000 


49000 


63000 


28000 


Decay time (ns) 


230 


25 


35 


40 


Peak emission 


415 


330-352 


358-385 


420 


Hygroscopic 


Yes 


Yes 


Yes 


NO 



5 Collection efficiency of the scintillation light 

photons plays another significant role in bar detector 
performance. Light collection performance is determined by 
the sensitivity and quantum efficiency of the photodetector . 
In the prior art, PMTs and PS-PMTs were used as photosensor 

10 devices. Recent innovations have improved the performance 
of photosensors. These developments have also provided new 
photosensors with superior characteristics. These include 
inorganic photodiodes (silicon drift, PIN, avalanche) and 
organic photodiode arrays. 

15 The conventional estimator of the Z position used in 

the prior art is as follows: 

z = k\^^\ + r (1) 

20 where k is a scale factor, r is an offset, El and E2 

are the total energy signals of the Dl and D2 PMTs, 
respectively. Since equation (1) is a linear estimation of 
the event position based on the non- linear ratio of the 
photosensor energy signals along the bar, inefficient and 

25 biased estimation is unavoidable. In addition, post 

processing of the bias correction technique is required. 

Positioning estimation is greatly improved by using a 
statistics based un-biased estimator for the invention. The 
novel estimator is follows: 
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z = arg min 



V 



J 



(2) 



where R is the ratio between the Dl and D2 energy 

E — E 

signal for a given event, i.e., R = — and jUr(z) and 

E x + E 2 

cr^(z)are the mean and variance of the ratio R at a given z 

location, respectively. The JUr(z) and must be pre- 

computed and stored in a look-up table. The ratio R is an 
index of the spatial resolution. 

Fig. 5 shows a block diagram of one example of a 
position decoding circuit, using a position decoding scheme 
with a crossed wire anode PS-PMT. The positions of the gamma 
ray interaction in the bar detector module in the X-Y plane 
are calculated based on the summed anode signals from both 
PS-PMTs mounted at the ends of the module. As shown, the 
individual signals from separate X and Y position anodes are 
summed and inputted to an Anger logic matrix or a registered 
network technique to calculate the centroid of the 
distributed electric current at both of the crossed anodes. 
The Z direction position is computed based on the total 
energy of the signals from the photosensors using the 
statistics-based positioning scheme explained in equation 
( 2 ) above . 

Comparison of Figs. 2(a) and 2(b) shows the improvement 
in spatial resolution using the proposed statistics-based 
positioning (SBP) algorithm over the conventional centroid 
estimator algorithm of equation (1) . Specifically, by using 
the SBP algorithm, the FWHM in a 5 cm long bar detector 
improved from 1.67 mm using the centroid estimator (Fig. 
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2(b)), to 0.91 mm (Fig. 2(a)). In addition, the SBP 
algorithm does not require any linearity correction. 

Another significant property in achieving an acceptable 
bar detector for nuclear medical imaging application is the 
5 surface property of the bar detector. The surface treatment 
of the scintillation bar detector influences three factors 
that determine detector performance : (1) light collection 
efficiency; (2) linearity of light collection as a function 
of position; and (3) the light collection ratio between the 

10 photosensors . Possible surface treatments may include the 

addition of reflectors, polishing (grounding) and etching of 
the long sides of the bar detectors. 

According to one preferred embodiment, the four long 
sides of the bar surface are grounded and sealed with a high 

15 reflectivity material such as Ti02 , MgO powder/paint or 

Teflon® tape/film for optical isolation and maximization of 
light collection performance of each scintillation bar. The 
surface condition versus expected spatial and energy 
resolution is shown in Figs. 6(a) and 6(b). Polished surface 

20 with reflector, polished surface with metal coating and 

grounded surface with reflector cases were tested with a 4x4 
mm 2 cross section by 5 cm long bar. In all cases, 
reflectivity of the reflector was chosen to be 98%. This 
simulation demonstrates that the grounded surface provides 

25 optimum attenuation of scintillation photons at the 

surfaces, maximizing the slope of ratio R, defined in 
equation (2), as a function of location Z. This yields an 
enhanced spatial resolution . 

In addition, the surface treatment provides a linear 

30 light response, which represents a uniform spatial 
resolution with respect to position. 
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As described above, the novel bar detector module -based 
gamma camera for medical imaging exhibits many advantages 
over the prior art, including: 

1) Three dimensional position sensitivity; 

2) Improved spatial resolution through use of statistics- 
based positioning algorithm, which also provides unbiased 
estimation of a gamma ray interaction position in the 
detector; 

3) Maximized position estimation accuracy and uniformity 
along the z -direction by using grounded surfaces and other 
surface treatments for the scintillation bar; 

4) Further enhanced resolution by use of new photosensors 
with improved sensitivity and quantum efficiency; 

5) Improved system performance through use of new high 
light output , fast decay scintillators ; and 

6) Ability to implement many diverse custom system 
configurations and associated performance and cost 
characteristics based on the selection of orientation and 
dimensions of the bars and modules, allowing for use in 
nuclear medicine imaging systems. 

The invention having been described, it will be 
apparent to those skilled in the art that the same may be 
varied in many ways without departing from the spirit and 
scope of the invention. Any and all such modifications are 
intended to be included within the scope of the following 
claims . 
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